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Abstract The production of the massive bosons Z® and 
could provide an excellent tool to study cold nuclear matter 
effects and the modifications of nuclear parton distribution 
functions (nPDFs) relative to parton distribution functions 
(PDFs) of a free proton in high energy nuclear reactions 
at the LFIC as well as in heavy-ion collisions (HlC) with 
much higher center-of mass energies available in the future 
colliders. In this paper we calculate the rapidity and trans¬ 
verse momentum distributions of the vector boson and their 
nuclear modification factors in p-tPb collisions at = 

63 TeV and in Pb-tPb collisions at -^snn = 39 TeV in the 
framework of perturbative QCD by utilizing three parametriza- 
tion sets of nPDFs: EPS09, DSSZ and nCTEQ. It is found 
that in heavy-ion collisions at such high colliding energies, 
both the rapidity distribution and the transverse momentum 
spectrum of vector bosons are considerably suppressed in 
wide kinematic regions with respect to p-tp reactions due to 
large nuclear shadowing effect. We demonstrate that in the 
massive vector boson productions processes with sea quarks 
in the initial-state may give more contributions than those 
with valence quarks in the initial-state, therefore in future 
heavy-ion collisions the isospin effect is less pronounced 
and the charge asymmetry of W boson will be reduced sig¬ 
nificantly as compared to that at the LHC. Large difference 
between results with nCTEQ and results with EPS09 and 
DSSZ is observed in nuclear modifications of both rapidity 
and pt distributions of and W in the future HIC. 
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1 Introduction 


The production of massive vector bosons in high- 

energy nuclear collisions has long been regarded as an excel¬ 
lent tool to probe the initial-state cold nuclear matter (CNM) 
effects in the ultra-relativistic heavy-ion collisions (HIC) at 
the LHC because Z^ jW^ do not participate in strong inter¬ 
actions and the mean-free-path of massive vector bosons (or 
more precisely, the leptons decayed from Z^ jW^ ) is rather 
long IMl. In proposed heavy-ion programs at the Lu- 
ture Circular Collider (LCC) at CERN ifl^ , and the Circular 
Electron Positron Collider with a subsequent Super proton- 
proton Collider (CEPC-SPPC) in China 1131, the vector gauge 
boson productions may still play a very important role in 
investigating the high-density QCD in the initial-state and 
making precise constraints on CNM effects and nuclear par- 
ton distribution functions (nPDEs) at very small x. We note 
that the massive vector bosons are produced very early at 
time ~ 1 /mz ly ^ 10^^ fm/c with a decay time smaller than 
0.1 fm/c, and a mean-free-path ^ 10 fm in the QGP at tem¬ 
perature 1 GeV II 14 II . Even though the hot/dense QCD mat¬ 
ter to be created in nucleus-nucleus collisions at the ECC 
and the CEPC-SPPC may have a much higher initial tem¬ 
perature and longer lifetimes Ifsll . the vector boson with the 
final state colorless di-lepton would still be nearly blind to 
the medium evolution 114, BS- 

Recent studies have shown that the heavy gauge boson 
production in nuclear reactions at the LHC energies could 
shed light on several CNM effects, especially the shadowing 
and anti-shadowing effects {for Z/W) and isospin effect (for 
111 , 3 - 23 ]. In the conceptual designs of ECC 1151 
and CEPC-SPPC ll24ll the center-of-mass energy for proton- 
proton collisions could reach up to 100 TeV, which may give 
the energy of about ■^snn = 63 TeV for proton-lead col¬ 
lisions and about = 39 TeV for lead-lead collisions. 
With much higher colliding energies available in future HIC, 
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Fig. 1 (Color online) The rapidity distribution of Z® boson in p+p 
collisions at yj = 7 TeV. The ATLAS combined data are taken from 
Ref. 1^ . 
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Fig. 2 (Color online) The normalized boson transverse momentum 
spectrum in p+p collisions at y/s = 7 TeV. The CMS combined data are 
taken from Ref. (S^l . 


the vector boson would be produced by initial partons with 
the much smaller momentum fraction x, in which region the 
shadowing effect will be more pronounced and sea quarks 
and gluons may play a more important role. It is of great in¬ 
terest to see how the massive gauge boson productions help 
us understand the CNM effects and impose stringent con¬ 
straints on nPDFs in future heavy-ion collisions. In this pa¬ 
per we study the nuclear modifications of the vector boson 
production in the future Pb+Pb collisions at = 39 TeV 
and p+Pb collisions at = 63 TeV. The numerical cal¬ 
culations are performed by using the perturbative quantum 
chromo-dynamics (pQCD) program DYNNLO ll2-5n incor¬ 
porating the parameterized nuclear parton distribution func¬ 
tions (nPDFs) sets EPS09 DSSZ HI and nCTEQ ||28, 


bution functions (PDEs) /(v,/r) and partonic cross section 
^ab^v+x^ii+x which could be computed with the perturba¬ 
tive expansion, 


Ozjy 


= Y. dXadXbfalA{Xa,P)fblB{xb,P-) 

a,b-' 

^ (^) - ]ab^V+X^ll+X, 


(I) 


This paper is organized as follows. In Section|2] we briefly 
introduce the framework of our calculation and then com¬ 
pute the nuclear modification ratios for boson rapidity 
distributions and charged lepton pseudo-rapidity dependence 
of IT^ production in Pb+Pb and p+Pb collisions in the future 
HIC. In Section[3 we discuss the nuclear modification of the 
transverse momentum distributions of Z** and bosons in 
the future HIC. We present our summary in SectionlH 


2 Nuclear modification ratio for vector boson rapidity 
distributions 

Before discussing the nuclear modification of vector boson 
production, we briefly introduce the theoretic framework of 
this study, which is the same as in reference in . In hadronic 
collisions, the differential cross section of vector boson pro¬ 
duction A + B ^ V +X ll +X through the Drell-Yan 


mechanism yC [jCould be calculated within the QCD-improved 
parton model ll3ll l32|] . In this approach, the cross section 
could be factorized into the convolution of the parton distri- 


where Xa(b) are the momentum fractions of the parton a{b) 
from the hadronA(B). As the baseline of the study, the hadron- 
hadron cross sections are numerically calculated with the 
program DYNNLO ll^ . which provides the ^(a*) and 
pQCD corrections to the leading-order (LO) cross section 
The ff{as) calculation includes real correc¬ 
tions and one-loop virtual corrections, and that at 
includes double-real corrections, real-virtual corrections, as 
well as two-loop virtual corrections lE^ . In addition, the 
parton distribution functions parametrization sets MSTW20- 
08 are used, and the renormalization and factorization 
scales are set at the boson mass niv- 

We plot in EiglUthe rapidity distribution, and in Eigl2] 
its transverse momentum spectrum, for proton-proton colli¬ 
sions at = 7 TeV. Eor rapidity dependence, both ff{as) 
and &{a}) pQCD calculations agree well with the ATLAS 
data Odll . which is measured in the invariant mass interval 
of the lepton pair as 66 GeV< mu < 116 GeV and also in 
the lepton transverse momentum region as pj > 20 GeV. 
The calculation at also give good description of the 

CMS transverse momentum distribution ll3^ in the region 
Pt > lO GeV. The CMS measurement is made in the fi¬ 
nal state phase volume defined by the invariant mass, trans¬ 
verse momentum, and pseudo-rapidity of the dilepton as: 
60GeV< m„ < l20GeV, p^j > 20 GeV, and Iq'] < 2.1. The 
DYNNLO program has also been shown to give good agree- 
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Fig. 3 (Color online) The nuclear modification ratio Raa of the Z** boson rapidity dependence for = 2.76 TeV (LHC, left panel) and 
ySiw = 39 TeV (future HIC, right panel) at the NLO. The LO results with EPS09 LO nPDFs are also plotted for comparison. 


ment with the CMS data on rapidity distribution, and with 
the transverse momentum spectra of and W measured by 
ATLAS m. 

To include several important CNM effects in Pb+Pb and 
p+Pb collisions, we phenomenologically utilize the parametriza- 
tion sets of nuclear parton distribution functions (nPDFs) 1126 - 
. In this work, we chose MSTW2008 proton PDFs 
and multiply them by the flavor and scale dependent factors 
Rf{x,ji) taken from EPS09, DSSZ, and nCTEQ to obtain 
the parton distribution of a bound proton in nu¬ 

cleus. With the isospin symmetry being assumed lEhl - E^ . 
the nuclear parton distribution of a bound neutron (.r, /i) 
could be obtained. In this paper only minimum-bias (MB) 
nuclear collisions are considered. 

Then we focus on the nuclear modihcation on vector 
boson production in the future HIC. Eirst we study that on 
Z® rapidity distribution in the PbH-Pb collisions at = 

39 TeV. The Z^ signal is dehned in the invariant mass win¬ 
dow of the hnal state lepton pair as 60 GeV <inu<\20 GeV, 
which is the same as the CMS experiments ISist]. The Z — y* 
interference is included in the calculation. We calculate at 
the order i^(oCj), the nuclear modihcation ratio RAA{y^) de¬ 
hned as 


RAAiy^) = 


da'^/dy^ 


nCTEQ nPDEs may suppress more than 30% of Z^ yield 
at the central rapidity, relative to proton-proton collisions. 
EPS09 and DSSZ give similar nuclear variations for the fu¬ 
ture HIC. The LO results with EPS09 are also plotted in 
Fig.El It is observed that the higher order corrections make 
negligible alteration on the 7?^^ at the LHC, but give a larger 
suppression on R^^ at future HIC. 

To give a simple picture of the difference between the 
two collisions energies, we resort to a lowest-order analy¬ 
sis d . In the LO pure electro-weak processes, the rapidity 
of Z'^ is related to the initial quark and anti-quark momen¬ 
tum fractions with the kinematic equation 


Xl7 = 




(3) 


{Ncoii)daPP/dy^ 


( 2 ) 


where the narrow width approximation (mu mz) is used. 
Equation (O shows that the momentum fraction of the ini¬ 
tial parton is inversely proportional to the colliding energy 
Eor example, one has x ^ 0.033 for mid-rapidity Z® 
production at 2.76 TeV, which falls into the anti-shadowing 
region of EPS09, whereas at 39 TeV x ^ 0.0023 is given by 
Eq. (O, that is small enough to enter the shadowing region 
of EPS09. The suppression in the 39 TeV Pb-tPh collisions 
is due to the shadowed smaller-x initial partons. To visually 
display the difference between the two energies, we plot in 
Fig.|4]the factory written as 


with (Ncnu) the number of binary nucleon-nucleon colli¬ 
sions 139]. 

In Pig. [3] both the results at the LHC and future ener¬ 
gies are plotted for comparison. Signihcant differences be¬ 
tween RAA{y^) at two energies can be seen. The LHC re¬ 
sults (left panel) with three nPDPs sets show weak nuclear 
modihcations (< 10%) in the studied rapidity region. En¬ 
hancements with EPS09 and nCTEQ at mid-rapidity can be 
observed, where suppression is given by DSSZ. However, at 
future HIC (right panel), strong suppressions can be found 
in the whole studied rapidity regime (|y| < 4.0). Especially, 


+ , (4) 

where Rf{x,jj.) is the flavor dependent nuclear variation fac¬ 
tor given by EPS09, DSSZ, or nCTEQ, and xi 2 is related to 
y^ by Eq. (|3. The subscripts v and s stand for the valence 
and sea quarks (no gluon at leading order), respectively. This 
R (y^, m) factor could reflect the nuclear modification on 
the LO partonic subprocess initiated with /v_i and fs, thus 
could shed light on the RAA{y^) fll . For the four processes 
plotted in Fig. |4l similar differences between the two colli¬ 
sion energies could be observed. Compared to those at the 
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Fig. 4 (Color online) The subprocess dependent factor 
for Pb+Pb collisions at both ~ 2.76 TeV and 39 TeV with fac¬ 
torization scale set at mz- 
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Fig. 5 (Color online) The contribution rates of various partonic sub¬ 
process as functions of in Pb+Pb collisions at both = 

2.76 TeV and 39 TeV calculated at the leading order. 


LHC, the PDFs in future Pb+Pb collisions are more depleted 
in the kinematic region of the mid-rapidity production, 
which results in a larger suppression of RAA{y^) in the future 
HIC than that at the LHC as observed in Fig. [2 We also no¬ 
tice that, although at the future HIC the My -|- Mj channel with 
nCETQ gives less suppression than those with EPS09 and 
DSSZ, their contribution combined with relatively large sup¬ 
pressions by other channels with nCTEQ, may give a more 
suppressed Raa^J^) than those with EPS09 and DSSZ. The 
underlying reason is that the Wy -|- Us channel gives relatively 
small contribution to yield in the future HIC, and we will 
discuss this issue in more detail in the following. 

We have mentioned that at LO Raa as a function of the 
vector boson rapidity mainly depends on the nuclear modi¬ 
fication on the quark and anti-quark distributions fill] , con¬ 
sidering the LO process qq^V. Actually anti-quarks only 
come from the nucleon sea, while quarks could be the va¬ 
lence or from the sea. To better understand the contributions 
to Raa given by some certain flavors, we calculate at LO 
the contribution ratio of partonic subprocess to the differ¬ 
ential cross section. Three primary kinds of processes for 
both LHC and future HIC are shown in Eig. 0 (l)My + «, 
and dv -I- dg', (2) Ug + Ug and dg -f dg and (3) i -I- i and c -I- c. 
Processes (2-3) completely depend on the sea quarks, while 
the process (1) contains valence quarks. Obvious difference 
could be seen between the two energies. At LHC energy, the 
processes with valence quarks in the initial-state dominate 
the Z® production. However in the future HIC their con¬ 
tribution will be reduced significantly, while the processes 
with pure-sea quarks in the initial-state become more im¬ 
portant. The results demonstrate that sea quarks play an im¬ 
portant role in the Z^ boson production, especially at the 


future HIC. A main reason is that, the sea quark densities 
increase very fast with the decreasing momentum fraction x, 
whereas the valence quark densities decrease sharply (more 
discussion could be seen in the [Appendix A| i. We show at 
the bottom of Pig. 0the corresponding momentum fractions 
of the incoming partons (xi: forward moving; X 2 '. backward 
moving). The strong suppression of Z'^ production in mid¬ 
rapidity region at the future HIC is lead by the shadowed 
nucleus sea quarks. It is noted that the nuclear modification 
of gluon density will also give contribution through high or¬ 
der corrections. 

We note that at the small-x region the uncertainties of 
nPDPs are rather large, thus the Raa given by different nPDPs 
may even merge with each other if error bar is included. 
That means the nuclear effect is very unclear at small-x, and 
the future Pb+Pb collisions will provide invaluable informa¬ 
tion to constrain the nPDPs for sea quarks and gluons in the 
small-x. Especially, the nuclear variations of s and c quarks 
densities, which has been constrained loosely by vector bo¬ 
son production at the LHC, should become important at the 
future HIC. 

The yield of Z^ boson in the future p+Pb collisions at 
y/sy/Ai = 63 TeV is also studied. Pigure0 shows the nuclear 
modification factor RpA as a function of the Z^ rapidity in 
the center of mass frame, predicted by EPS09, DSSZ and 
nCTEQ nPDPs. Results for LHC energy (g/sjm = 5.02 TeV) 
are also shown for comparison. Since the nuclear effects 
come from the colliding lead nucleus, the asymmetric rapid¬ 
ity dependence of RpA could be observed at both LHC (left 
panel) and future HIC (right panel) energies. Por the LHC 
energy, enhancement can be seen in the backward region. 
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Fig. 6 (Color online) The nuclear modification factor RpA {y^) for p+Pb collisions at ^/smn ~ 5.02 TeV (LHC, left panel) and ^/snn = 63 TeV (fu¬ 
ture HIC, right panel) at the order 


while suppression is given in the forward region (related 
to anti-shadowing and shadowing, respectively llill ). For 
future energy, except the enhancements in very backward 
regime (—4.0 <y^< -3 .0), persistent suppressions are given 
by EPS09, DSSZ, as well as nCTEQ. 


The LO analysis is simple for pH-Pb collisions, and one 
could study the RpA from the nPDFs factors by 

replacing the momentum fraction X 2 (related to the nuclear 
parton) with y^, according to Eq. (O. In Fig. |7] we plot the 
Rfiy^jjJ.) factors of certain flavors for both LHC and future 
HIC. Compared to the LHC energy, the future pH-Pb colli¬ 
sions generate more suppression in the rapidity region > 
—2.0 due to the shadowed initial parton from the lead nu¬ 
cleus with smaller-x. Though for s quark density in the LHC 
backward region, and Wy quark density in forward region of 
both energies, the nuclear modifications with nCTEQ show 


distinct behaviour, their contributions to rapidity distri¬ 
bution will be rather small in the corresponding region. To 
see that clearly, three kinds of contribution of the partonic 
subprocess initiated by certain nuclear flavors are illustrated 
in Fig. 0 It is observed that for both LHC and future HIC, 
the valence quarks’ (uy and dy) contribution decreases with 
the rapidity (smaller-x), but the sea quarks’ increases. Com¬ 
parison between the two energies implies that the nucleus 
valence quarks densities give negligible contributions with 
the increasing energy (smaller-x). The suppression in the fu¬ 
ture pH-Pb collisions is to a great extent the results of the 
shadowing effect on the sea quarks. The contribution of the 
s and c quarks will also become non-trivial (more than 30% 
in the mid-rapidity at LO). The Raa and RpA as functions 
of the rapidity would become a good probe of the nu- 
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Fig. 7 (Color online) The flavor dependent nPDFs factors Rf (jZ, jl ) 
for p-rPb collisions at both =5.02 TeV and 63 TeV with factor¬ 
ization scale set at mz. 


Fig. 8 (Color online) The LO rapidity dependent contribution rates 
for partonic subprocess initiated with various nuclear flavors in p-rPb 
collisions at both y/s^/M = 5.02 TeV and 63 TeV. 
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Fig. 9 (Color online) The nuclear modification factor RaaChO f™' boson production in Pb+Pb collisions at y'S/wv = 2.76 TeV (LHC, top panel) 
and y/Sfltl = 39 TeV (future HIC, bottom panel) at the order 


clear modification on small-x sea quarks in future heavy-ion 
collisions. 

For completeness, we also study the nuclear modifica¬ 
tion on the W boson production. Usually, the charged lepton 
pseudo-rapidity, instead of the boson rapidity, is measured 
for W boson production in experiment with the neutrino 
in final state (W Iv). We calculate the Raa and RpA as 
functions of the charged lepton pseudo-rapidity for W pro¬ 
duction in Pb-fPb collisions at = 39 TeV and p-fPb 
collisions at ^/snn = 63 TeV, respectively. The transverse 
momentum redon of the charged lepton is chosen to be 
p'j > 25 GeV HB]- 

For nucleus-nucleus collisions (Fig. IS, asymmetry be¬ 
tween Raa{W^) and due to the nuclear isospin 

effect, can be observed. To be specific, because of the ex¬ 
istence of neutrons in nuclei, the nuclear reaction may en¬ 
hance the production of and reduce the yield of 1T+, 
relative to those in the proton-proton collisions. From the 
panels (d) and (h), one can see that the separation between 
1T+ and is large at the LFIC, and becomes small at the 
future Flic, especially in the mid-rapidity region. Actually, 
isospin effect is related to the parton distribution asymme¬ 
try (m(x) ^d{x), and m(x) ^ d{x )), which originates mainly 
from the asymmetric ujd valence quark distributions in a 
nucleon {(uud) for proton, and {udd) for neutron). As is dis¬ 
cussed before, at the future energy the initial parton tends to 
come from the smaller-x regime, where the sea quarks domi¬ 
nate the valence quarks. Thus, at smaller-x the isospin effect 


becomes weaker (more detail can be found in Appendix A[ ), 
and nuclear modifications on VT+ and become more 
symmetric. It is also noted that, the nuclear effect of the total 
W i=W^ -\-W^) production, is very similar as that of the Z** 
rapidity dependence at the same collision energy. The differ¬ 
ences among the Raa with three nPDFs sets are also consis¬ 
tent, for various particles (IT, W^, W^, and Z^) at the same 
collision energy. Three nPDFs sets predict similar nuclear 
modifications at the LFIC (panels (a-c)), whereas nCTEQ 
gives much stronger suppressions in the future FIIC (pan¬ 
els (e-g)). Since isospin effect is very weak at the central 
rapidity, the suppression on the W boson production at the 
future HIC is mainly due to the shadowing effect on the nu¬ 
clear sea quarks. 


For proton-nucleus collisions (Fig. doll, the nuclear vari¬ 
ation on is again the result of the isospin effect versus 
(anti-)shadowing effect. At the LHC, one can observe that 
the isospin effect results in a separation between 1T+ and 
in the backward direction (larger x of the nucleus), and 
shadowing effect suppresses their production in the forward 
region (smaller x of the nucleus). However, the isospin ef¬ 
fect is rather weak at the future energy, and suppressions on 
both VT+ and in a wide range can be observed. 


The charge asymmetry observable £/ of IT^ production, 
defined as 


Alw+ -Nw- 
Nw+ +A^w- ’ 


( 5 ) 
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Fig. 10 (Color online) The nuclear modification factor Rp^{ri^) for W boson production in p+Pb collisions at = 5.02 TeV (LHC, top panel) 
and y/SivAf = 63 TeV (future HIC, bottom panel) at the order ^(a.,). 


is also studied for p+p, p+Pb and Pb+Pb collisions at both 
y/sim = 5.02 TeV and = 39 TeV, and are plotted in 
Fig. [TT]as a function of the charged lepton pseudo-rapidity. 
One can observe the charge asymmetries at the LHC show 
quite distinct behaviors in three different colliding systems, 
whereas the three curves at the future energy lie near the 
horizon and are close to each other. It has been found that the 
nuclear variation of the charge asymmetry observable £/ of 
is mainly due to the isospin effect, and not sensitive to 
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Fig. 11 (Color online) The charge asymmetry as a function of the 
charged lepton pseudo-rapidity in p+p, p+Pb and Pb+Pb collisions for 
both the LHC (yTjvV = 5.02 TeV) and future HIC (yTj^ = 39 TeV). 


the other nuclear effects fill] . Therefore it is understandable 
that the isospin effect will become weak in the future HIC, 
as shown in Fig. [TT] 


3 Nuclear modification ratio for vector boson transverse 
momentnm distribntions 

In this section we investigate the CNM effects on vector bo¬ 
son transverse momentum spectrum, which may be used to 
get access to quite different kinematic regions from that by 
the boson rapidity production iflll] . In this paper we mainly 
focus on the massive boson production with > 10 GeV, 
where DYNNLO provides an excellent baseline description 
for pj distribution in elementary hadron-hadron collisions. 

First, we study the nuclear modification on boson 
transverse momentum distribution in future Pb+Pb collisions 
at ^/snn = 39 TeV. Besides the same Z® invariant mass inter¬ 
val as used in the last section, the rapidity range \y^\ < 2.0 is 
also chosen. The factor 7?^, as a function of boson trans¬ 
verse momentum pj, is calculated at the order ^(cti). Al¬ 
though the differential cross section will be changed a little 
by higher order corrections, Raa defined as the ratio of the 
cross section of A+A to that of p+p would not be very sen¬ 
sitive to the higher order corrections or the variation of hard 
scales Predictions by EPS09, DSSZ, and nCTEQ, 

for both LHC and future HIC, are shown in Fig. [12] Ob¬ 
vious difference between the two collision energies can be 
seen. Appreciable enhancements are given by three nPDFs 

























Fig. 12 (Color online) The nuclear modification ratio of the transverse momentum distribution in Pb+Pb collisions at y/s^ 
2.76 TeV (LHC, left panel) and = 39 TeV (future HIC, right panel) at the the order ^(a^). 
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Fig. 13 (Color online) The flavor dependent nPDFs factors Rf(pj,yi) 
for Pb+Pb collisions at both .Js^ = 2.76 TeV and 39 TeV with the 
factorization scale fixed at mz. 


Fig. 14 (Color online) The transverse momentum dependent con¬ 
tribution rates of various partonic subprocess in Pb+Pb collisions at 
both -y/Siw = 2.76 TeV and 39 TeV calculated at the LO. 


sets in the region 10 GeV< pj < 250 GeV at the LHC, 
whereas considerable suppressions are given by them in the 
region 10 GeV< pf < 150 GeV for the future Pb+Pb colli¬ 
sions. The nCTEQ nPDFs supports stronger nuclear modifi¬ 
cations than the other two in the studied pr regime. 

To well understand the differences among the nPDFs 
sets and also between the two collision energies, one can 
again perform the analysis at LO. At this order, we obtain 
the kinematic relation 

PT + Jpl+ml 

___y_ /z'. 


at mid-rapidity with the narrow width approximation. Then 
according to Eq. ©, we replace the initial parton momen¬ 
tum fraction x with pj and plot the flavor dependent factor 


Rf{pj,p) of nPDFs in the Fig. [T3] For both LHC and fu¬ 
ture HIC, those factors of EPS09, DSSZ, and nCTEQ are 
shown. Compared to the nuclear effect at the LHC, more 
suppression on the parton distribution could be seen in the 
region pj < 200 GeV for the future nucleus-nucleus colli¬ 
sions, except for the nCTEQ u valence quark distribution. 
Because the related initial parton momentum fractions are 
different (smaller at future energy), shadowing effect dom¬ 
inates the future Pb+Pb collisions whereas anti-shadowing 
effect is shown up at the LHC. 

To see clearly how the nuclear modihcation on each fla¬ 
vor comes into the observable RAA{pf)’ we again calculate 
the the contribution rates of the various partonic subpro¬ 
cesses for both the LHC and future HIC. Three kinds of 
contribution are plotted in Fig. [141 (1) Partonic subprocesses 
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Fig. 15 (Color online) The nuclear modification factor of the ZP boson transverse momentum spectrum for p+Pb collisions at 
5.02 TeV (LHC, left panel) and y^sjvjv = 63 TeV (future HIC, right panel) at the order 


with valence quarks in the initial-state, such as My + itj or 
dv + g', (2) Pure u,d sea initiated subprocesses as Us + Us 
and ds -I- d/, and (3) Gluon initiated partonic processes with 
u,d sea like u^ + g. One could observe that processes of 
Type (1) give the dominant contributions at the LHC en¬ 
ergy, but their contributions decline significantly in the fu¬ 
ture HIC. On the other hand, processes of Type (2) and Type 
(3) and other sea quark initiated processes dominate in the 
future energy collisions. In the previous work dii on 
the nuclear effects on pj spectra of gauge bosons at LHC, it 
has been shown that the gluon’s contribution is predominant. 
Here we demonstrate that the impact of the nuclear modih- 
cation of the sea quark distribution, on the Raa{Pt)’ com¬ 
parable to that of the gluon, in the future nucleus-nucleus 
collisions. For example, the large enhancement in the re¬ 
gion 100 GeV< pf < 300 GeV at the LHC with nCTEQ, 
is mainly due to the gluon nuclear modification but not the s 
quark (also enhanced a lot by nCTEQ), while at future HIC 
the shadowing of both gluon density and sea quark (u,d,s,c) 
distributions in nCTEQ give signihcant suppressions to pr 
distributions of Z®. 



Fig. 16 (Color online) The LO partonic subprocess contribution rates 
as functions of for different nuclear flavors participating in p+Pb 
collisions at both V5 aw = 5.02 TeV and 63 TeV. 


The nuclear modification factor of the boson trans¬ 
verse momentum distribution in the future p-fPb collisions 
at s/si^ = 63 TeV is also calculated at i^(oCs), as shown 
in Pig. [15] Similar to that in PbH-Pb collisions, the suppres¬ 
sion due to the shadowing effect could be observed at the 
pH-Pb in the future HIC. The nuclear modification in p-fPb 
collisions is weaker than that in Pb-fPb collisions because 
in pH-Pb only one colliding object is lead nucleus. The par¬ 
tonic subprocess contribution rates are calculated for both 
LHC and future p-fPb reactions. Pig.[T6|shows the rates for 
three kinds of process, each of which is initiated by certain 
flavors in the lead nucleus. One can see that for both of the 
two energies, processes initiated with sea quarks or gluons 
give the main contribution. Especially at future energy the 


valence quarks’ contribution is marginal, the nuclear effects 
are brought in by the sea quarks and gluon. 

Next we turn to study the nuclear modification on the W 
boson transverse momentum spectrum. The W boson pro¬ 
duction is calculated in the final state phase space defined 
by the center of mass frame pseudorapidity and transverse 
momentum of the charged lepton as jTj^l < 2.5 and pj < 
25 GeV 101. In Pig. [T7|we plot (Pr ) PbH-Pb collisions 
at both = 2.76 TeV and ^/snn = 39 TeV at the or¬ 
der ^(ctj). The curves for and are separated from 
each other by the isospin effect, between which lies Raa for 
total W (= + W^). The separation increases with the 

W transverse momentum (increasing with x), and decreases 
with collision energy (decreasing with x). Although pro¬ 
duction obviously depends on the isospin effect, that of total 
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Fig. 17 (Color online) The nuclear modification ratio Raa(Pj) for VT+, W and total W{= VT+ + W ) productions in the Pb+Pb collisions at 
= 2.76 TeV (LHC. top panel) and = 39 TeV (future HIC, bottom panel) at the order 


W receives a Z'^-like nuclear variation mainly because of the 
(anti-)shadowing effect. One can observe the Raa of bo¬ 
son are very different for two collision energies. It is because 
that the isospin effect strongly enhances the production 
at the LHC, while the shadowing effect is more pronounced 
than the isospin effect in the region <100 GeV at the 
future HIC. It should also be noted that, the suppression on 
at the LHC is due to the isospin effect, but that at the 
future energy is caused by the shadowing effect. 

The nuclear modification ratio RpAiPr ) ^Iso calcu¬ 
lated for the future pH-Pb collisions at y/s^N = 63 TeV, and 
the results are shown in Fig. [18] Results for the LHC energy 
are also shown for comparison. Same trends as in PbH-Pb 
collisions could be seen in pH-Pb. It should be noted that al¬ 
though both isospin effect and (anti-)shadowing effect are 
weaker in pH-Pb as compared to PbH-Pb, the shadowing ef¬ 
fect still surpasses the isospin effect in the <100 GeV 
region and thus suppresses the yield in the future HIC. 

Next we study the charge asymmetry observable s/ as 
a function of the W transverse momentum. The results for 
pH-p, pH-Pb and PbH-Pb collisions at both y/s^N = 5.02 TeV 
and y/SNN = 39 TeV are plotted in Fig. ED One can observe 
that the nuclear modification of £/ due to isospin effect in¬ 
creases with pp and decreases with the collision energy. 
Actually .e/’s nuclear modification also increases with its 
value in proton-proton collisions which is related to the par- 
ton distribution asymmetry ratio r„j(T) and rpj{x) (defini¬ 


tion can be found in Appendix Aj l. Figure. [19] demonstrates 
the isospin effect will be reduced in the future HIC. 

From the numerical results we note that the differences 
among various nowadays nPDFs sets are considerable. As 
an example, we calculate the ratio of the Raa{Pt) the 
LHC (y/SNN = 2.76 TeV) to that at future energy (y/SNN = 
39 TeV) for total W production, as 

R^^^ip^ + Ap^) 




R^^'-{p^,-Ap^y 


(7) 


The results for p^^ = 90 GeV are shown in Fig. ]20] One 
can observe that EPS09 and DSSZ predict similar ratios of 
Raa at the LHC and future HIC, with very flat Ap^ depen¬ 
dence. However, nCTEQ gives a very different result: the ra¬ 
tio is obviously larger than those with the other two nPDEs, 
and increases faster with Ap^. The measurement on nuclear 
modification of the vector boson production in the future 
heavy-ion collisions will provide very useful constraints on 
nPDEs. 


4 Summary 

We investigate in this paper the initial cold nuclear mat¬ 
ter effects on the heavy vector boson production in the fu¬ 
ture nuclear collisions. Within the framework of perturbative 
QCD, we numerically calculate the nuclear modification ra¬ 
tios of the boson rapidity distribution, the charged lep¬ 
ton pseudo-rapidity dependence for W production, as well 
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Fig. 18 (Color online) The nuclear modification factor RpA as a function of the W boson transverse momentum for the p+Pb collisions at y/s^N - 
5.02 TeV (LHC, top panel) and.yijvAf = 63 TeV (future HIC, bottom panel) at the order 



Fig. 19 (Color online) The charge asymmetry as a function of the W 
boson transverse momentum in p+p, p+Pb and Pb+Pb collisions for 
both the LHC (y/iAw = 5.02 TeV) and future HIC (y(siw = 39 TeV). 



Fig. 20 (Color online) The ratio ^(Apf ) of the Raa {Pj ) at 
2.76 TeV to that at y/sjUl = 39 TeV for total W production. 


as jW^ transverse momentum spectra. Both Pb+Pb col¬ 
lisions at y/SMN = 39 TeV and p+Pb collisions at = 
63 TeV are studied, where the CNM effects are included by 
utilizing the parameterized nuclear parton distribution func¬ 
tions sets (EPS09, DSSZ and nCTEQ). 

Different nuclear modifications from those at the LHC 
energies could be found. Eirst, the nuclear shadowing ef¬ 
fect would play a more important role in massive gauge bo¬ 
son production in the future HIC. As consequences, both 
the (pseudo-)rapidity dependence and transverse momen¬ 
tum distribution of the boson would receive suppressions 


in a wide regime. Second, the nuclear isospin effect would 
be rather weak in the future nuclear reactions, which could 
be clearly observed in the production. Eor example, the 
production in the future HIC is even suppressed in the 
<100 GeV region, though at the LHC it is enhanced by 
isospin effect. 

These very distinct nuclear modifications at the future 
HIC with those at the LHC result mainly from the smaller 
momentum fraction x carried by the initial partons at the 
fixed rapidity or pj. At the LHC, the momentum fraction 
of the initial-state partons may fall in the anti-shadowing re- 






























12 


gion (for EPS09), however in the future HIC it may be seen 
in the kinematic region of the shadowing effect. Further¬ 
more, the calculations of the contribution rates for various 
partonic subprocesses show that the valence quarks gradu¬ 
ally lose their relevance with the increasing colliding energy 
because of their sharply depleted densities in the smaller-x 
region. With the highly suppressed u and d valence quark 
distributions, the valence quark density asymmetry is mi¬ 
nor, and the isospin effects would become much weaker in 
the future HIC. 

In addition, large differences among results with sev¬ 
eral parametrization sets of nPDFs could also be seen at 
the future HIC. For example, results with nCTFQ predict 
more suppression for the vector boson yields than that with 
FPS09 and DSSZ in future HIC. The theoretical understand¬ 
ing of nuclear modifications of PDFs in nuclear collisions is 
far from satisfactory, and the massive gauge boson produc¬ 
tion in the future HIC could provide an excellent probe of 
the shadowing effect on the nuclear sea quark distribution 
as well as gluon density by imposing new and precise con¬ 
straints on the nPDFs at rather small-x. 


Appendix A: c(s{x) and r„_d(x)/r^ j(x) 


] 

0.8 

0.6 

0.4 

0.2 

0 

0.2 

0.1 

0 

- 0.1 


'■-(a) / 

- /'■ 

: / \ 

/ 

/ V 

(b) ' / J 

^^s(X) - 

- valence^' \ : 

- sea / - 

,(c} / 

JudiX) 

(d) ; 

rudix) : 

' .. .'ifl0. ih®. 

. . . 11^ .. 1 .. 


io-‘ 


10-^ 

X 


1 10 -^ 


10-^ 

X 


2.76 TeV 
•■• 5.02 TeV 
■39- TeV. 
.63 TeV 


Fig. 21 (Color online) Panels (a,b): the c{,s(r) of valence and sea for 
u and d quarks in a proton. Panels (c,d): the parton distribution asym¬ 
metry ratio r„ ;/(x) and r- j{x). The factorization scale is set at the Z® 
boson mass rnz. The marks at x-axis are corresponding to the different 
colliding energies. 


C? J(■^)) should be emphasized. Actually, the parton distri¬ 
bution asymmetry is to a large extent caused by the valence 
quark (mv, dy) density asymmetry, as a result the tenuous va¬ 
lence quark density would give a much weaker r^d (or 
in the small-x regime. 


We calculate the c{,^s{x) of the valence and sea quarks in a 
proton with the MSTW2008 PDFs Q as 


Cf{x) 


C{ix) 


Mx) 

fvi.x)+fs{x) 

fs{x) 

fv{x)+fs{x) 


l-Cif(x), 


(A.l) 


where f = u,d. The results of u and d quarks are plotted in 
the panels (a) and (b) of Fig. |2T] respectively. One can see 

f 

that for both u and d, valence quark Cy (x) increase with the 
parton momentum fraction x, and decrease with the collision 
energy. 

The parton distribution asymmetry r^dix) and r^j{x) 
defined as 


rud{x) 

rudi^) 


u(x) —d(x) 
u(x) -\-d{x)'’ 
u{x) — d(x) 
u(x) -h d{x) 


(A.2) 


are also plotted in the panels (c) and (d) of Fig. |2T] One 
could observe that the parton distribution asymmetries are 
small in the low-x region. As a result, the isospin effect 
would be moderate in that region. We also mark at the x- 
axis the points corresponding to different energies according 
to the relation x ^ my / y/s^N- The marked values in pan¬ 
els (a) and (b) are obtain with my = niz and those in (c) and 
(d) with niy = my/. The correlation between the c({x) of 
valence quark and the parton asymmetry ratio r„ j/(x) (and 
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